Frame-dragging and the kinematics of Galactic- Center stars 
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ABSTRACT 

We calculate the effects of frame dragging on the Galactic-Center stars. As- 
suming the stars are only slightly relativistic, we derive an approximation to the 
Kerr metric, which turns out to be a weak field Schwarzschild metric plus a frame 
dragging term. By numerically integrating the resulting geodesic equations, we 
compute the effect on keplerian elements and the kinematics. We find that the 
kinematic effect at pericenter passage is proportional to (a(l — e^))^^. For known 
Galactic-center stars it is of order 10 m/s. If observed this would provide a 
measurement of the spin of the black hole. 

Subject headings: gravitation — relativity — stellar dynamics — Galaxy: nucleus 



1. Introduction 

The center of the Milky way is a very interesting region. It contains a massive black 
hole (MBH) of ~ 3 x IO^Mq. The central parsec contains thousands of stars . For a small 
but a growing number of these, due to the relatively close proximity of the MBH and short 



orbital periods, the orbital parameters h ave been accurately measured (ISchodel et al.l 12003 
Chez et al.ll2005l : lEisenhauer et al.ll2005l ) . So me of the star s have pericenter velocities as 



high as a few percent of c. Hence, as shown in IZucker et al.l (120061 ) the general relativistic 
effect of 0(/3^) should be observable. 

At 0{l3^) general relativity predicts a new effect, which is that a spinning black hole 
drags the surrounding space-time along with it. Under the rotational frame-dragging effect 
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(also known as Lense-Thirring effect), the frame of reference with minimal time dilation is 
one which is rotating around the object as viewed by a distant observer. If this effect could 
be observed for GC stars, then in principle the spin of the MBH can be measured. A method 
based on the orbital dynamics of the GC stars is more direct than the usual approach which 



requi res modeling the effect of spin on the accretion disk (see for example iNarayan et al. 
2008h . 



The effect of O{0 ^) t erms on the kepleria r i elern e nts and on astrometry were discussed 
bv ljaroszvnskil JlQQsh and lPrarile fc Mathewsl J200oh . Iwilll J2008h goes on to consider 0{[3^) 
as well. The resulting astrometric effects are so small that they can only be observed on 
stars which are closer to the MBH than the observed GC stars. 

In this paper we concentrate on the effects of the 0{(3^) terms on the kinematics of 
the GC stars. Traditionally the effect of relativistic perturbatio ns on orbital dy namics have 
been studied either u sing post-newtonian cel estial mechanics ( IWeinbergI [1972) or pseudo- 
newtonian equations (ISemerak fc Karaslll999l ). We adopt a different and conceptually sim- 
pler approach. We do a low velocity perturbative expansion of the Kerr metric and then 
numerically integrate the resulting geodesic equations. 



2. The Model 

Our starting p oint is the Kerr metric in Boyer-Lindquist co-ordinates (see for example 



Misner et al. 197i 



ds"" = (dt -ssinHd^Y + ^ ((r2 + s'')d(P - sdi? + ^dr^ + p'de'' (1) 
pZ \ J pi \ 'IS. 

where 

A = - 2Mr + ^ + cos^ Q (2) 

Here s is the spin of the black hole and M is the mass and G = c = \. 

The Kerr metric describes the space-time outside a rotating black hole. The metric itself 
is complicated and is difficult to solve even numerically. So we make some approximations 
to make the solution easier. We consider the case where 

~ 1/r (3) 

This is equivalent to saying that the system is not very relativistic. To agree with Eq. ([3]) 
we assume r is 0(e~^) and (ir, rdQ and rdc^) are all 0(e). We then replace 

r — > e~^r (4) 
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dr edr (5) 
de e^de (6) 



e 



3. 



(7) 



in the Kerr metric, put M = 1 and Taylor expand up to O(e^). We get the following metric: 
ds'' = - (l dt^ + + e^dr^ + e\^de^ + eV^ sin^ Odcj)^ - ^ sin^ edtd<P (8) 



which is equivalent to a weak Schwarzschild field plus a frame-dragging effect. The above 
metric is only valid for systems with /3 <^ 1. In particular, it is not valid for null geodesies 
as light does not satisfy Eq. Q. 

Applying the Euler-Lagrange equations, we get the following geodesic equations for t, 
r, 6 and 0. 

1 ]i+—ri + 0{e') = Q (9) 

2e2\ . . {2 p 2se^ 

1 + — f = + r sin^ Ocj? - — + e^— + sin^ Oicj) (10) 

r^e + 2rre = — sin2et0+ (11) 

2se^ /sin^ 6^ • \ 

sin^ ^0 + sin 2000 + 2r sin^ 0r0 + ri - sin 200t - sin^ 0i = (12) 

r \ r / 

The leading order terms give the Newtonian equations, terms up to O(e^) gives Schwarzschild 
equations and terms up to O(e^) gives the frame- dragging effect. We remark that e as used 
here is just a label for keeping track of orders. For numerical work, we set e = 1. 



3. Results 

We now proceed with the numerical integration of the geodesic equations fl9]- [T2|) . For 
simplicity we choose the unit of length to be the gravitational radius GM/c^ of the central 
MBH (~ 5 X 10^ km) and velocity to be in units of speed of light, which makes the unit of 
time to be around 17s. 

We choose a, e, J, VL and uj as the free parameters, where a is the semi-major axis of 
the ellipse, e the eccentricity, / the inclination with respect to the spin axis of the black 
hole, VL the longitude of ascending node and uo is the argument of perihelion. We start the 
integration at apocenter and integrate from t = to t = 27ra^/^ which is exactly equal to one 
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newtonian orbital period. We then re-compute a, e, /, fi, at the end of one integration and 
hence find the effect on the keplerian elements. Relevant formulas for the initial conditions 
are given in the Appendix. 

We actually integrate three different sets of geodesic equations: 



• The equations for the newtonian case (includes only zero order e terms). The orbital 
elements are completely unchanged. 

• The Schwarzschild case (terms up to e^). There is a pericenter shift in this case which 
we call Aus- 

• The frame dragging case (terms up to e'^). There is a shift in the node AQ^ and a 
further shift AcUfd- There is no net effect in a, e and /. 



We have also computed the velocity differences between these cases. 

We examine the numerical data to verify the known parameter dependencies of the 
various relativistic effects or find an empirical formula. All the effects turn out to depend 
upon 

p = a{l-e'^) (13) 
which is the square of the angular momentum in gravitational units. 
First the well known expression for the pericenter shift 

Qtt 

Acus = — (14) 
P 

has been verified numerically as shown in Fig. [TJ 

Next we consider the frame-dragging precession, for which we empirically infer: 

An,, = ^ (15) 

and 

Aufd = -127TS COS (16) 

as Fig [2] verifies. For a plane perpendicular to the spin axis of the black hole (cos/ = 1), 
Afifd + Aa;fd cor responds to the total pericenter precession which we see is equal to — Svrp^^/^ 
(cf. Section 3 of [Weinberg et al.ll2005l ). The physical shift of the apocentre is ~ a{l + e)AQfd- 
As viewed from a distance of 8 kpc this translates into an astrometric shift of 



Aa = a(l + e)Afifd x 4yuas 



(17) 
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It should be noted that we have calculated the change in keplerian elements after one 
complete revolution around the MBH. Dividing by the orbital period 27ra^/^ gives the mean 
rate of chaiige. Do ing so in Eqs. [HI [TS] and Uni gives expressions matching Eqs. 6, 4 and 5 of 
Jaroszynskil (jl998l ). 



We now consider kinematic effects. Fig.[3]shows the comparison between the Schwarzschild 
and frame- dragging effects on the three velocity components. The effect is a fraction of a 
kilometer per second during pericenter passage for typical GC star orbit parameters. Note 
that the velocity components Vr = r, vq = r9, v^f, = r sin 6*0 are actually derivatives with 
respect to the conformal parameter. Finally, we consider the maximum velocity difference 
between the frame-dragging and Schwarzschild effects (A^d) and between Schwarzchild and 
Newtonian effects (AV^). Fig. H] shows the empirical relations: 



AK 



8e 



p 



,3/2 



AVf 



fd 



—8Aes cos J 
p2 



(18) 
(19) 



All these empirical relations are fairly accurate for parameters typical of GC stars but 
not exact. On the figures the points do not exactly lie on the line which suggests that these 
are leading order effects. 

Table 1 gives the values of the various relativistic effects according to our empirical 
formulas for a sample of GC stars and two binary p ulsars. For the GC star s we estimated p 
using Eq. f[T5]) . from the a and e values tabulated in lEisenhauer et al. For t he binary 

pulsars we derived p using Eq. (IT^ and the precession rate given in IWilll (120061 ) . We see 
from the table that the GC stars are more relativistic than binary pulsars. The advantage 
is that the binary pulsars have very short orbital periods (less than a day) and hence we get 
many more orbits. 



4. Conclusions 



We see that the maximum kinematic effect in known GC stars AVfd is of the order of a 
few lO's of m/s during the few weeks around the pericenter passage. Although this level of 
accuracy is difficult to achieve for GC stars, it is not implausible. Extrasolar planet searches 
regularly reach an accuracy better than 1 m/s (ILovis et al.ll2006l ) and new technol ogies for 
radial velocity measurements may be able to obtain precision as high as 1 cm/s (ILi et al. 
2008h . 



There are also two theoretical problems which remain to be solved. 
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First, an accurate calculation of the redshift as a function of time is required (as it is 
the observable quantity), rather than velocity as a function of time as calculated here. Both 
the kinematic and gravitational redshifts are involved. We do not know of any approximate 
method for calculating the redshift in this case, as our approximate metric is not valid for 
light. It may be necessary to calculate null geodesies in the full Kerr metric. 

Second, the relativistic effects have to be separated from the newtonian effects of other 
masses, such as nearby stars, gas and dark-matter clouds. These could overshadow the 
frame-dragging contribution, especially since some newtonian d ynamical processes i n the 



GC region can be un expectedly strong because of resonances (iGiirkan &: HopmanI 12007 



Lockmann et al.ll2008l ). However, newtonian perturbations from other masses would not 
give the distinctive time dependence in the kinematics that frame-dragging does (Fig. [3]). 
Hence, we can be optimistic about dis-ent angling frame-dragging from all the newtonian 
effects. 
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A. Evaluating the orbital elements 

The numerical integrations in this paper are done in Boyer-Lindquist coordinates, whereas 
the resuhs are presented in terms of Keplerian orbital elements. To convert between them, 
we use standard relations from celestial mechanics. In practice, we only need to use the 
classical formulas at or near the apocenter, so we will treat r, 9, (p as ordinary spherical polar 
coordinates. 

For the initial conditions, we need to set up a star at apocenter with given a, e, /, Q, u. 
Since the relativistic effects are minimal here we set i = 1. We start by defining a temporary 
cartesian coordinate system, centered at the black hole, but oriented such that the star is 
on the X axis with velocity along +y. In other words, we put the star at 
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where 



r. 



apo 



a(l + e) ^;apo = J —— (A2) 



1 + e 

Applying the rotation 

R,{Q) R^{I) R,{uj + tt) (A3) 

(right operator first) gives the position and velocity in the reference cartesian system. We 
then convert to spherical polar coordinates. 

For the inverse process at the end of an integration, we start by computing position r 
and velocity v in cartesian coordinates. We then compute the specific angular momentum 
h and the Runge-Lenz vector e. 

r 

h = rxv e = vxh (A4) 

r 

The inclination / and the longitude of the ascending node Vt are simply a way of specifying 
the orbital plane, and we have 

Q = arctan(/ij^, /ij;) + ^ I — arctan ^-y//i^~+^, (A5) 

Here arctan means the two argument form also called ataii2. For the argument of the 
perihelion, we consider the Runge-Lenz vector (which is a vector having magnitude e and 
pointing towards the pericenter) in the orbital plane 

ecostc; ' 

esintu I ^ R,{-I) R^{-9) e (A6) 
and the left hand side gives uj and e. 
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Table 1. Relativistic effects in GC stars and Binary Pulsars 



Object 


P 






s cos I 
(radiansj 




(km/s) 








AVs 
(km/s) 


Astrometric shift 


SI 


8.6x10* 


5.0x10" 


7 


-1.5x10" 


6 


-1.2x10" 


4 


2.2x10" 


4 


0.034 


0.27 


S2 


6.8x10^ 


2.2x10" 


5 


-6.7x10" 


5 


-4.7x10" 


2 


2.8x10" 


3 


3.7 


4.9 


S8 


1.1x10* 


1.1x10" 


5 


-3.2x10" 


5 


-1.9x10" 


2 


1.7x10" 


3 


1.9 


6.5 


S12 


1.3x10* 


8.7x10" 


6 


-2.6x10" 


5 


-1.4x10" 


2 


1.5x10" 


3 


1.5 


4.5 


S13 


4.4x10* 


1.4x10" 


6 


-4.1x10" 


6 


-5.1x10" 


4 


4.3x10" 


4 


0.10 


0.40 


S14 


6.4x10^ 


2.5x10" 


5 


-7.4x10" 


5 


-5.8x10" 


2 


3.0x10" 


3 


4.4 


10. 


B1913+16 


2.9x10^ 


8.1x10" 


8 


-2.4x10" 


7 


-1.9x10" 


5 


6.5x10" 


5 


9.5x10"^ 




PSR J07370 


2.3x10^ 


1.2x10" 


7 


-3.5x10" 


7 


-4.3x10" 


6 


8.3x10" 


5 


1.9x10"^ 





3 

<] 




Fig. 1. — Pericenter precession Acj^ due to O(e^) terms, plotted against for a sample of 
orbits with different a, e, /. The line shows the classical result of Eq. f[T^ . 
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Fig. 2. — Frame- dragging precession Aflf^ and AtUfd due to 0{e^) terms,plotted against the 
dependencies ( fTSi) and ( fT6l) respectively, for a range of orbits as in Fig. [T] 
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Fig. 3. — Comparison between the Schwarzschild (dashed curves) and frame-dragging (sohd 
curves) contributions to the three velocity components f^, Vq and v^. The time starts from 
the apocenter. For this orbit a = 3 x 10"^, e = 0.90 and / = 0.5. 
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Fig. 4. — Peak kinematic contributions AVg and AVfd, plotted against the dependencies f|T8l) 
and (fT9l) respectively. 



